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Ie INTRODUCTION 


During 1946-1949 sevoral articles have been published by the 
University of Texas giving half-hourly observations of both tempera- 
sure and wind at various levels from the ground up to and including 
268 fect for periods of aporoximately one daye This suggested tho 
problem of determining the coefficient of eddy conductivity using the 
observed lapse rate data; and, the coefficient of viscosity by using 
the observed wind profile data.s These two coefficients were found, in 
1936, by Sverdrup D2] to be equal in tho atmosphere, but other eddy 
diffusion coefficients wero found by him to be unequale On the other 
hand, Petterssen end Swinbank [oj 3 in 1947, suggest that the former 
coefficient is more nearly related to the latter by the factor 65 in 
the free atmospheree The question appears to require another critical 
examination; and, with the availability of new ond fairly accurate data, 
it seemed to be a problem well worthy of study. In actuality, the problen 
of determining the diurnal variation of the coefficient of eddy conductivity, 
K, near the ground proved to be so complex that this was all that time per-= 
mitted in this dissertation. It will be necessary to answer the question 
of the equality of the two coefficients at some later date. However, the 
deternination of values of the coefficient of eddy conductivity and its 
diurnal variation, determined purely from heat transfer considerations, 
is of some practical interest as a check on previously determined valuese 
For example, Stewart [i] indicates a possible range of 10? em seeml to 
10? om@ seo" for K within the surface layer, depending upon the stability 


of the layers 
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This study gives explicit values of K, with respect to local civil 
time, corresponding to a particular synoptic situation, vnat Lor the 
period 27-28 September 1948 in the vicinity of Manor, Texas. Another 
important result of this study proved to be the determination of ratos 
of radiative heating and cooling of the layers of air near the ground at 
various times covering a full day. This yielded some surprisingly large 
valuese These results could be used to solve a problem proposed by 
Brunt fod concerning the magnitude of the coefficient, Kes of radiative 
diffusivity. Brunt originally computed Kp to bo about 650 em? soc7l but 
in 1940, using the new absorption coefficients of Dennison, Ginsberg and 
Veber [4] he concluded that this value of Kp probably should bo increased 
thirty-fold. The water vapor radiative transfer problem of this thesis 
used the Elsasser Atmospheric Radiation Chart which, in turn, is based 
primarily on the new absorption coefficients mentioned abovee Actually 
Kp was not computed in this study but the results of this study form a 
basis which would permit the calculation fairly sinplye 

It was found also that carbon dioxide plays a surprisingly large part 
in the radiative heat transfor processes near the ground where rather large 
lapse rates are found. In the preparation of his radiation chart, Elsasser 
made the assumption that the carbon dioxide contribution to the not flux 
could be neglectcd in comparison to that due to water vapore It turned 
out, in this study that the radiative contribution to the rate of change 
of temperature due to carbon dioxide was in the neighborhood 1/3 to 1/4 
that due to water vapors Elsasser's remark, concorning the effect of carbon 
dioxide, "that even thin layers radiate as blaci: bodies within these bands" 


was showm to be definitely invalid. In fact, using a formla for carbon 
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dioxide absorptivity due to Caliendar [3] » the computed absorptivities 
ranged anywhere from 2% to 20% for layers near the earth's surface, de- 
pending upon their thicknesse 

Surprisingly large values of turbulent temperature-change por half 
hour were obtained, in this study, at three inches above the grounde for 
example, there was a turbulent rate of cooling of 10°C per half-hour at 
1300 local time. It was found, however, that even such large values wore 
negligible insofar as affecting the magnitude of the coefficient of eddy 
conductivity in the surface layer. In other words, whenever tho lapse 
rates were significantly different from dry adiabatic, it was found that 
the turbulent temperature-change term of the partial differential equation 
of eddy conduction was negligible in comparison to the turbulent heat flux 
through either boundary of a slab of air. Thus the surface layer extended 
at least to 80 feet throughout the period of this investigation and during 
the night extended to about 200 feet. 

Other results deduced in this study are the following: The coefficient 
of eddy conductivity in the surface layer appears to follow a power law of 
the form K = K, 2", with m71 the exact value of m depending in some way 
upon tho stability. loreover in this layer, K takes on its maximm value 


at about 1300 LCT and its minimum value at about 1900 LCT. 
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IIe THE NATURE OF THE DATA 


le The Data. 
Data utilized in this research was taken from Report Noe 29 com= 
piled by the Electrical Engineering Research Laboratory of the University 
of Texas [6] e This report contains recorded measurements of temperature 
and vapor pressure for selected levels from the ground up to and including 
288 feet. The soundings are recorded at half-hourly intervals throughout 
the 24 hour period chosen to represent a particular symoptic situatione 
Measurements of wind speed and direction are recorded for selected levels 
up to 307 feet. -The sounding site is on level ground under cultivatione 
‘The authors of [6] give the following statements of accuracy of the 
data. 
{a) Tho temporatures are accurate to within 1°C over tho 0-50°C ranse. 
(b) The d=ninute averaging period used to odtain tomperatures is not 
sufficient to remove all turbulent fluctuations; however, they do not deviate 
more than *,2°C from the 30-minute mean temperature This turbulent fluctu- 
ation of the 6-minute mean as compared with the 30-minute mean will be called 
“che gust error". 
(c) The Aerovanes show a uniform response with an accuracy of +2%. 
ror a detailed description of the data and the field installation, tho 
reader is referred to the report of the data [6y . 
2e Selection of Data for this Study. 
The data for the period September 27-28, 1948 wes solected for this 
study. A high pressure cell dominated the area and gave conditions of 


horizontal homogenietye 
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5, 185, 235 and 288 feet, 
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III. THEORETICAL AND EMPIRICAL CONCEPTS 


le The Vater Vapor Radiative Transfer Problem. 

Elsasser [5] gives a method for determining the net flux due to 
the water vapor content of a layer of atmosphere through the use of his 
radiation charte The corrected optical depth of water vapor, u,; ina 


column of air of unit cross-section is computed by using the formla: 
oo yf , > 
U = ~ Jooo > t 7000 L\e ea (361) 


When only temperature and moisture data are available at given levels, 
prossures for these levels may bo computed from the surface pressure by 


uso of the integrated form of the hydrostatic equation: 


gaz-k, i (An. p, - tn pe). (3.2) 


For the evaluation of net flux into or out of a slab of air, utilizing 
areas on Elsasser's Atmospheric Radiation Chart, Elsasser recommends the 
measurement of these aroas by means of a plenimeter, The difference of the 
net fluxes for a slab is then given in gram-calories per three hours by this 
method and is converted to rate of change of temperature of the slab by means 


of the equation: 


£2 F, = -¢,(aT) Ae : ee 


(6) 
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In this equation, Fp is the net upward flux at the top of tho slab, while 
F, is the net upward flux at the bottom of the slab. For further details 
concerning this method, reference should be made to Hlsasser'’s study, Heat 
Transfer by Infrared Radiation in the Atmosphere [5] e 
2. The Carbon Dioxide Radiative Transfer Probleme 

Callendar [3] gives formlas for the absorptivity in terms of the 
quantity of carbon dioxide for various wave bands of a parallel beam of 
radiation. His formmlas for the equivalent thiclmess, c, of carbon dioxide, 
and for the monochromatic absorptivity of carbon dioxide, in a columm of air 
of thiclmess Z, pressure Po (in atmospheres) and temperature T (degrees 


absolute) are: 


m 
u 


8.7Z pa T ; (3.4) 


AH 1-1/[itncrpy]. (3.5) 


In equations (304) and (305), ¢ is the length of absorbing path in ems of 
carbon dioxide at I].f.Pe and is based on air having a normal carbon dioxide 
content (the partial pressure of carbon dioxide in the eir now being .0032 
atmospheres). Tho constants n and x are given for different limits of the 
intense absorption band centered at 15 microns. 

Blsasser [5] states that, in the case of square-root absorption, a 


Slab of thiclmess u is mathematically equivalent for diffuse radiation to 


(7) 


a linear columm of length 1.78 u for beamed radiations Tho condition of 
square-root absorption is satisfied for slabs of air of thiclmess 4 to 100 
meters according to Panofsky [10] e For thinner layers, the factor 1.78 
should be increased slightly; but the error introduced for such layers will 
not be great if the factor 1.78 is usode Inasmuch as diffuse radiation 
occurs in the atmosphere, slab thiclmesses mst be increased by the factor 
1.783 which amounts to an increase in the equivalent thiclmess, c, of 
equation (304) by the factor 1.78. Utilizing this corrected equivalent 
thickness in the absorption formula (305) allows this formula to be applied 
to slab thicimesses of the atmospheres 

In this study, the rate of change of temperature of layers of air near 
the ground due to carbon dioxide as well as water vapor is to be determined.e 
Furthermore, since the Atmospheric Radiation Chart of Elsasser is to be used 
for the rate of temperature=change resulting from water vapor, consideration 
mst be given to his treatment of carbon dioxidee LElsasser [5] considers 
that the absorption by carbon dioxide occurs, for the most part, in the range 
752 cml to 584 om71, that is, between 13.3 and 17el micronse He further 
considers that "very thin layers radiate as black bodies" and, therefore, 
neglects the radiative effects of carbon dioxide in comparison to those of 
water vapor, for problems concerning net fluxe The chart, however, does 
give a separate representation of the black body flux contained within the 
strong band interval 1303 - 1721 micronse Callendar [3] gives the values 
N= 03) and X m 055 as suitable for determining the absorptivity by (5) 
within a band of limits 15.0 - 17.0 microns. It is assumed that the values 
of monochromatic absorptivity, A, determined from (34) and (3.5) are appro- 


priate for the carbon dioxide limits of Elsasser, 1363 - 17.1 microns. 
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Knowing the ground temperature and the mean temperature of the 
layers ek tmosphere under consideration, the flux through the vop 
boundary of each layer can be computed. For exemple, for the first slab 
above the ground, the equation of radiative transfer, sometimes called 


Schwarzchild's equation, gives 
Eom (1-A,) Exo tA, Baa » (306) 


where [I is the flux through the top boundery of the layer, A, the mono-= 
chromatic absorptivity of the layer, Yo che black body flux in the band 
1365 - 17.1 microns at the temperature, T, of the ground (or the flux into 
the lower boundary ) and Ey, is the black body flux at the mean temperature 
T of the layer. Strictly speaking Schvarzchild's equation (3.6) is true 
for monochromatic intensity, or beamed flux, of radiation. Actually 
Elsasser [5] integrates the Sclwarzchild equation in order to obtain the 
flux emerging from an arbitrary slabe He finds that ae integration in= 
volves the computation of integrals E; i, (x) = Sy y-* en dy where 
Pzsec Yy; values of which are given in cole However, ho also shows that 
these computations imply the increase of tho absorbing path by the factor 178 
discussed abovee It is therefore assumed that Senvwarzchild's equation holds 
for fluxes E, Eyo, end Eqy With the value of absorptivity A corresponding to 
the absorbing path c replaced by 1079 c, a result which Elsasser justifiese 
For slabs above the first, equation (3.6) also permits the computation 
of upward flux E through the top boundary. However , for such computations 
Exo must be taken to represent tho flux into the bottom of the slab, which 
was thet out of the top of the slab immediately below, EQL would still 
represent the black body flux, in the band 13.3 = 17.1 microns, corresponding 


to the mean temperature of the slab under consideration. 
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In order to apply the same procedure to the problom of computing 
fluxes for the dowmmard radiation of carbon dioxide as for the upward 
radiation, it is necessary to know the flux entering the top boundary of 
the highest layer. The top boundary is 268 feet for this study. Through- 
out the period used in this study, whenever a ground inversion oxists, this 
height is very near the top of the inversion. The slope of the sounding 
above the sround inversion during this period is quito similar to that of 
the standard atmospheree This remark is, on the whole, true also of the 
slope of the sounding above the 288 feet level during the hours in which 
no ground inversion is present. Therefore, the atmosphere from 10 kilo- 
meters is divided into slabs of approximatoly 100 meters thicilmess and 
having the appropriate temperature according to the standard atmospheric 
lapse rate. In treating equation (36), slab 1 is to be considered the 
topmost slab, slab 2 that immediately below, etc. It will be shown below 
that the amount of flux from above 10 kilometers produces negligible 
effect at the 288 foot level. Let flux E, enter the top of slab 1. The 
fluxes leaving the bottoms of the next several slabs are given by the 


following table. 


Slab Flux through the lower boundary 
1 1 ay ay ; 
2 ml Ej +tma+b, (307) 
4 nn Ey +nma #nb tec. 


In (3-7), 1, m, n..... ete. correspond to tho factors. (1-A)), (1-4, ); 
(1-Az); eecee OtC., respectively, of equation (3.6); a, b, ¢ correspond to 


the values Ay By), Ay Byos Ay Bas eves Gtoe OF Goucion (3.6). 
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Since, using equations (3.4) and (3.5), the average value of 1, m, 
Nl, escoes toe is of the order 0.5 for a layer of thickmess 100 meters; 
and B, is of the order of 10 cal om@? (3 hr.)"*, it turns out that the 


contribution of the initial flux E, is negligible in treating equation 


° 
(3.6) for the dowmward directed flux. ie may therefore arbitrarily set 

=, = 0 in determining the flux through the bottom of the slabs and still 
obtain the correct result at 288 feet. 

Finally, since we now may determine the amounts of upward and dowmiord 
directed flux at every level, one can compute net-upward, net-downward and 
consequently net-outward directed flux, Fo - Fy for every consecutive pair 
of levels. This determines the radiative heating, or cooling, of the slab 
due to carbon dioxide radiation, and the actual rate of temperature-change 
can be computed, as before, from equation Ga) 

3 Heat Transfer by IMolecular Conductivity. 


The upward flux of heat due to molecular conduction is given by ‘he 


formila: 


oT 
eS al (3.8) 
where @is the coofficient of molecular conductivity of air, and has the 
value &@ = 549(10)7/ cals en™+ sec™? 77+ at 0°C, but is a function of tem- 
perature (c.f. for oxample, pe 44 of [2] )e 
It is necessary now only to compute upward fluxes Fg and Fal for the 
various slabs for which lapse rates ere know, and again use equation (3.3) 


to convert net flux out of a slab into a corresponding rate of change of 


(11) 


temperature. In actual practice it is normally found, since “% is vory 
small, that in order to register a measurable net flux, the changes in 
lapse rates along the vertical mst be quite extreme , and this can be 
expected only within the lowest two or three slabse 

4. Rate of Change of Temperature by Eddy Conductivity. 


The formula for the observed rate of temperature change is: 


ot = (=). a oT) + (22). +(P). (3.9) 


oT 
where (2 is due to water vapor, 


¥). is due to carbon dioxide, 
| due to molecular conduction, 
ot - duo to eddy conductione 


The observed rate of change of temperature per halfhour is given for 
each layer by the data used in this studye As a result of methods previously 
described, the only unlmovm of formla (38) is the rate of change of tempera~ 
ture due to eddy conduction. 
5e Tho Coefficient of Eddy Conductione 

The neteoutvard flux for a layer of atmosphere due to turbulent con- 
duction may be expressed as: 

ke ES Zz) - Ie (52) = =¢ (gr) & 17 fa (3010) 
iAP C ZL 2 fra P prot a 
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n (3010) Ky and Ky are the coefficients of eddy conductivity at the top 


and bottom of any slab of atmosphere, respectivelye Making the substitu- 


tions - oi a a Prn 


= Be GE G8 


both approximations being valid to within 1% for the layers between the 


surface and 288 feet, we obtain: 


2.93 eles + = Mesh ry) ‘ iit WH) sa 


In (3e11), all items labeled with subscript 1 refer to measurements made 
at the center of slab 1, all items with subscript 2 refer to the center of 
slab 2; p, is the mean of p, and Poe Thus (3e11) determines the rate of 
change of temperature approximately at the boundary of slab 1 and slab 2. 
In general, slab 1 will denote any lower slab and slab 2 the slab immed- 
iately above slab 1. Thus (3011) may be used as a recursion formla to 
determine all values of K; after K,_ 1 has been determined. The only limi- 
tation on its use arises when -25 - x » that is, when the lapse rate 
of the slab is dry adiabatic, or, not significantly difforent from x e 

In using (3.11) as a reeursion formmla, it is necessary to determine 
& value of K at some specific level to permit determination of all other liye 


In order to do this without making too many assumptions, a value of K at 


(13) 


some standard level can be determined by reference to the theory of 
nonentum=transfer within the surface layer. The coefficient of eddy 
viscosity can be determined at a standard level, say 16 feet, and it can 
be assumed to give a representative value of the coefficient of eddy con= 
ductivity. The particular values of K will then depend to some extent on 
monentun-transfer theory, but the comparative values reflect only the 
measurements of eddy conductivity. Recent results of Lottau [9] are 
dravm upon for the purpose of finding a value of K at some standard level. 


‘According to the Prandtl theory of momentum-transfer*, onc obtains: 


k= (24). (3.12) 


Ou 


In (3612) 1 is the moan mixing length and OZ is the observed shear of 


v 


the horizontal wind, both measured at a given level. Lettau [9] gives the 
following formulas which permit the calculation of the coefficient of eddy 


viscosity within the surface layer: 


t= f/ +33/ 4] Ve =k (Z+Z,) , (3413) 


where 


OZ ° (314) 


*Haurwitz, B. Dynamic Meteorology. Equation (742). McGraw-Hill 1941. 
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In (3613) and (3014), tho subscript “a” indicates an isotropic or adiabatic 
' 8 
surface layers Tis the mean temperature of the slab for which 42. hes 


Ua 


been determined. The value of may be determined by the well-known 


logarithmic law 
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In (3015) ond (3215), all variables are understood to refer to the surface 
layere For example, 

T, is the constant value of the stress, 

U, is the horizontal wind speed, 

1, is the mixing length, 

k is the Von Karman constant, k = 45, 


Z. is the roughness parametere 


(@) 


OU, 
The values of Z, and \2 can be obtained from a plot of oe when 
adiabatic conditions are founde Observations of Lettau [3] indicate that 
such conditions exist at sunrise and sunset. in average u, profile, there- 


fore, can be dravm based on both sunrise and sunset data in order to detcr- 


mine Zo and the mean value of eS 


for the daye 


(15) 
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IV. RATES OF CHANGES OF TEMPERATURE 


le Rate of Change of Temperature Due to ‘viater Vapors 

Radiation transfer due to water vapor was computed for each layer of 
the soundings chosen from the data by the method described in Chapter III, 
Section 1, of this study. This method requires a Imowledsge of the specific 
humidity at each level. The data fo] gives the vapor pressure at cach 
level. These vapor pressures were converted to specific humidity by use 
of the Kiefer Multi-Prossure Hysrometric Chart. Since equation (3.1) -re- 
quires the pressuro-thiclmess, Ap, of all layers it was necessary to com= 
pute the pressure at every level by equation (302). To do this, the sur- 
face pressure must be kmowm. The surface pressure was computed from the 
sea level pressure recorded for Austin, Texas on tne surface map at the 
middle of the period under study, 27-28 September 1948. This was done by 
reversing the procedure of the reduction of pressures to sea level, knowing 
the height of the base of the tower (525 feet MSL). 

A (u,2) relationship was computed for all levols of the soundings and 
plotted on the Elsasser Chart. The difference of the net fluxes of each 
slab of the soundings was obtained by making the required area measurements 
With the aid of a planimeter. Figures 1 and 2 show schematically the areas 
that represent the differences of net flux, for one slab of atmosphere for 
tho following cases: 

(a) The slab lies within a ground inversion (Figure 1). 

(b) The temperature decreases with height (Figure 2). 

In Figures 1 and 2, the reference levels L, and L, denote respectivoly the 


2 
bottom and top of the slab, both taken to lie below tho 286 foot level. 
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Consider next Figure 1. ‘The (u,2) curve extending above Ly mst 
eventually cross thet above Lo because the former curve eventually corres- 
ponds to a sreater integrated optical depth of water vapor. Suppose this 
"cross-over" point is labeled Ae At 2200 LCT, the complete radiosonde was 
available to supplement the 2200 micrometeorological sounding. It was found 
that the cross-over points for the (u,‘) plots above the slab levels Ly and 
Lo occurred from about 300 feet for the lower, thinner slabs to about 800 
fect for the higher, thicker slabs. loreover, it was found that the area 
computation made by terminating the (u,T) curves above Ly and Lo at points 
B, and By corresponding to the 288 feot level was in error only by 2-73 
when compared with tho areas given by the complete (u,T) curvese This small 
error is due to a tendency of positive and negative areas below the isotherm 
ByBy to cancel each other. 

Bone iaes next Figure 2, representing typical daytime conditionse No 
exact comparison is known here of tho error made in neglecting the area 
below and to the right of the isotherm ByB5- Io cross-over point occurs in 
this case so there is no tendency for cancellation of positive end negative 
arease For the lower, thinner slabs little error could occur because the 
(u,T) curvos above the slab boundaries are in peneral nearly coincident at 
By and Bre The greatest error should occur in the slab 185-235 feet and the 
error should be progressively smaller in the slabs immediately below. The 
neglocted area, if appreciable, should cause a reduced amount of cooling or 
increased warming in the midday hours in the slab 185-235 feet. Actually 
Table 1 shows no evidence of reduced coolins of the slab 185-235 feete For 
tho ratio of the 1300 LCT cooling of this slab to those of the threo slabs 
imiediavely below is almost identical with the corresponding slab ratios of 


the 2200 LCT sounding, whose values arc knovm to be correcte 
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(Times given in LCT.) 
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Three~hourly Rate of Change of Temperature Due to ‘iater Vapor 
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Table 1 sives the values obtained for heat transfer due to water 
vapor in cal. (3 nr. )74, Table 2 gives the three=hourly rate of chance 
of temperature for each layer of each sounding, computed from Table l by 
use of equation (303). 
2e Rate of Change of Temperature Due to Carbon Dioxide. 

The differences in upward-directed flux due to the radiative transfer 
of carbon dioxide, in the band 1303 ~ 17-1 microns,were computed for each 
layer by the use of equation (366). The differences in dowmvard-directed 
flux wereootained from the method of Chapter III, Section 2. 

The computations of dovnnvard-directed flux took into account that the 
mid=period sounding was, on an average, ten degrees warmer throughout than 
that of the standard atmosphere, although the slope was assumed to be that 
of the standard atmosphere. This amounted to increasing the average emission 
from the layers by the factor 10%, which factor is the percentage increase 
in flux Dy, of equation (3.6) due to the increase in the moan temperature, 
tm» Of each layer compared to the standard atmospheric temperature. The 
adjusted value obtained for the flux through the top boundary of the layer 
286-235 fect was 18.79 cal em™@ (3 nre)72, This valuo was used as a constant 
value, for the period of this study and computations of flux differences for 
the layers of each sounding under study were made using equation (366)- The 
values of black body flux of emission for the carbon dioxide band 13.3 = 17el 
microns, at the mean temperaturo of oach slab were taken from the top of the 
Atmospheric Radiation Chart. The values of monochromatic absorptivity wore 


computed using equation (304) and (35). 
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results of Table 3 and equation (3.3 
_ = = . 


The tenperature=chenge values of Tables2 and 4 were addec alrebraically 
a ae oe” 

yielding a resultant radiative temperaturo-change for the slabse Since 

these temperature-changes are to be compared with the observed temperature-~ 

@ 

change at the various boundaries it was necessary to adjust them to these 

boundaries. This wes done by considering the radiative temperature-changes 

as existing at the midpoints of the slabs and then interpolating to the 

slab boundaries on log-log paper. These results are presented in Table 5, 


wherein all three-hourly rates have been reduced to half-hourly rates»by 


dividins by 6. 
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ed by the use of 
le 
The given in the observed data, for each 
+ 
< was considered as passing through the 
-_ = 
center of each layer. Consequently, the temperature-chanre computed 


oe 


é 
from the value of the flux differencos for two consecutive layers must 


be considored as bein at the boundary between the two layers, This 

resulted in a rate of temperature-change at the boundary of each slab 

without the necessity of interpolations The results obtained for rate 

of temperature-change due to molecular conduction are siven in Table 6. 
(Times siven in LCT) 


Boundary 02CO0 0500 0800 1100 1300 1600 1900 2200 


Gafeet) @e) (9c) (cc) (cc) () @&) 8) 


05 =090 -1490 -1660 7420 8.0 = 4410 057 = -2 260 
15 -l2 -216 209 BOT 045 wey ec] 0.0 
320 206 207 205 -.06 


Half-hourly Rate of Temperature=change Duo to 
Molecular Conduction 


TABLE 6. 


de Rate of Chance of Temperature Due to Eddy Conduction. 
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Those results follow immodiately by the use of the observed 9x 
and tho values of Table 5 and 6 in equation (3.9). The values for the 


rate of terperature-change due to eddy conduction are given in Table 7. 
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Ve RESULTS AND CONCLUSIONS 


le Pethod of Computation of the Coefficient of Eddy Conductivity. 

As indicated in Chanter III, Section 5, it was decided to devermine 
for all soundings evaluated a first value of K from the thoory of turbu- 
lent transfer of momentum, For this purpose it was necessary to make use 
of observed winds in order to evaluate Fa for use in eouation (3212). 
Two winds were siven consistently in the data which were always within the 
surface layer, one at 12 feet and one at 41 feet. For each of the eisht 
soundines studied, these tio values of the wind were plotted arainst log Z 
on semi-log paper, and a straisht line joining the two points was dravme 
The value of Va in the vicinity of 16 feet was then read off as being 
the value Ay corresponding to a height increment of 100 centimeters 
at 16 feet. The reasons for choosing the 16 foot level as the standard 
reference level are two-fold: 

(a) The level had to be within the lowest levels for which wind 
data are available. 

(b) The level had to be within the lowest slab consistent with (a) 
for which lapse rate data is available. This is necessary since, in 
equation (3013), the mixing lencth, 1, involves ea a (24 + x) : 

In satisfying roouirement (b), it was necessary to use the slab 12-20 
feet, which was considered to define oz at 16. feet. 

In several cases however, it was necessary to shift to the level 255 
feet. This necessity arose when an obviously inconsistent lapse rate was 


indicated in the layer 12-20 feet. This shifting to the level 23.5 feet 


(29) 


amounted to dealing with the over-all lapse rate in the layer 12-35 feet, 
which then geve , Ta at 2355 feets ee 
aie 

In order to evaluate the mixing length, 1 (the remaining rm of 
equation (3.12)) at the standard level, a kmowledce of the corresponding 


value 1,, the mixing length in an isotropic atmosphere, is necessarye 


a? 
According to Lettau ro] isotropic or adiabatic turbulence is nost nearly 
realized around sunset or sunrise, Thus, in order to derive the adiabatic 
velocity=profile for the 24-hour period under investigation, an average of 
three half-hourly wind-speed reports at both sunrise and sunset was obtained. 
Then, plottin« mean wind-speed arainst lor Z on semi-los paper for the two 
values of Z a straight line was obtained, Ficure 5, which according to 
Lettau should represent the wind profile for Manor, Texas under adiabatic 
conditions for 27-28 September 1948. The straight line obtained apparently 
does represent a close approximation to the true adiabatic atmosphere, since 
Tigure 3 indicates a roughness parameter Z, = 2 ems, which checks rather well 
with values quoted by Stewart [2] for fallow fields. MKnotting the value of 
the roushress parameter nermits the computation of le from equation (3013). 
The value of Sia at 16 feot was then read directly off the velocity 
profile of Firure 3. "Sines 16 feet 2 490 ems this can be done quite simply 
by reading off the wind-speed differences between 600 and 400 ems and divid- 
ins by 2C0 cms. This calculation leads to 415 sec~4, a value which may be 
checked nathematically using the theoretical equation of the’ straight line 


of Figure 356 This is done belows 
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The equation of the wind profile line of Pigure 3 can be written 
as in (3.15).Using the mean value of U, at the 12 foct level (366 ems), 
kom 045, 2, = 2 cms, Ug = 306 mps, we obtain ae" = 32.8 om sec7+, 


Then from equation (16), 


1 
2) no * (32.8) (Gas-) = «15 secon}. 
Ua 


Similarly at 23.5 feet the value of WA was found to be: 


Ou » os a! Rel 
WA = rag (5208 la = ike) sec ° 


At these two levels 1, = 220 and 323 coms, respectively. Values of K at 
the reference level may now be determined once the values of OZ > oe 
u 
® "Ta are Imowne ‘These results ere tabulated below in Table 8 for the 


various times of day: 


* o* * 
0200 0500 1100 1300 1600 1900 2200 
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Computation Values for Determination of K at Reference Level 


TAREE 6. 
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OO values of K were not computed since 


anaes 


ner non-representative cr dry adiabatic* lapse 
rates so that the heat transport equations (3.11) yield no solutions. 

Yor all soundings, other than the O8CO soundinr, it was possible to 
obtain solutions of ees: (3e11) using the values of K of Table 8 as 
a starting voint. The values of K fer 0200; 0560, 1100, 1900, 1600, 1900 
and 2200 LCT at the indicated elevations are listed in Table 9, and shorm 
graphically in Plate I. Im Table 9 the occurrence of blank spaces indi- 
cates the existence of a non-ropresentative lanse rate in the slab having 
the indicated height as its midpoint. It was therefore necessary to con= 
solidate adjacent slabs, obtaining a value of IK corresponding to the new 
midpoint, as discussed earlier in Section 1 of thie chapter. 

*Actually, isovnendiea lapse rates were not significantly cifferent from 


thé dry adiabatic for many of the lavers under study, when the gust error 
of ¥ .2°C is considered. 
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se 6 ad 
1900_ 2200 


K(10)3___—-K(10)_—sK (10)? 


026 004 036 0099 
307 066 716 
2ge5 ss 300 ea 1604 my 41 2.3 
405 a5 a4 5 150.6 073 ad 
: (eS aie 
| 940 469 509 20. 4o6 leo Oe] 
16.0 64.4 306 ere) 16. 
2365 13. 28. 41. 
| e155 500. 63.6 5 el Wed 
| eee) = 11s 236 460. 1700. 356 £0 Ves 
cn i 28. 300. 130. 32. 9.0 
§2.5 120. 
95 20 Be 2 le 136 20% 
7.5 10. 183 3506 25 5508 
woeO 0 lee 206 28. 840. 
0.0 17. 62. 2106 9206; 
oe 6 «(18 87. 190. 1500. 


The Coefficient of Dddy Conduction, K, (em sec™*) 


TABLE 9. 


(34) 


¢ - 
In order to = moro cloarly the acti of K with t 


and time, graphs of Bektecs elevations were 


variation of K at + 


Feet 


025 


at 


Diurnal Variation of "kK" 


4. 


igure 


I 


(36) 


o 
c 
1s 
.<] 
oa 


Feet 


me" Bec 


bh 


ion o 


Teese 


. 


1 Ver 


iurna 


D 


Pigune 5. 


(37) 


at 30 Feet 


Kk" 


mor ! 


iatio 


Diurnal Var 


Figure 6. 


(38) 


ion of "K" at 60 Feet 


e}euts 


1 Var 


aturna 


D 


Figure 7. 


(39) 


at 125 Feet 


gon of "7a 


t 


jurnal Varia 


D 


8. 


igure 


ig 


(40) 


iscussion of Xesnlts, 


t- 


2 
Lettau [9] has defined the surface layer essentially as the lazer 
in which the turdulent timo-rate of change of a property is noslisiole 


in comparison with its turbulent flux. In the solution of equations (3011) 


~ 


for the various levels, it is evident that the value of the lapse rates be-= 


come exceedingly critical as the value of the lapse rate approaches the dry 
adiabatic u, - For, if the maximum cust error 4,2 °C occurred in opposite 
senses at adjacent levels, it would cause an error in observation of the 
lapse rate by sevoral tines y » dependine upon the layer thickness. Thus 
it follows that the nost representative values of K may be expected where 
vhe lapse rate XY is largest compared to x and especially where it is of 
the order of 50-100 (cr more) times x « In order to emphasize this de# 
pendence upon lapse rate, let us consider the values of the lanse rate in 
the various layers for the 0200 and 1300 LCT sovndinrs. ‘These are showm 
in Table 10. In this table, the abbreviation "ner." indicates a non-~ 
representative lapse rate; therefore, necessitating a consolidation of 
the two adjacert layers. 

Notice, from Table 10, that at 1400 LCT the lapse rates were strongly 
super-adiabatic near the ground but not simnificantly different from the 
dry adiabatic above 80 fect. At 0200 LCT the lapso rates were more nearly 
uniform than at 1300 LCT. A comparison shows that the 1300 LCT lapse rates 
were numerically sreater than those for 0200 LOT between the surface and 6 
feet, but, above 6 feet, the early mornins lanse rates in seneral exceed 
those of mid-day. "or nurposes of comparison, the value of the dry adiabatic. 


lapse rate is x, = 095 x 1074 °C per come 
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In general, it was found that when ¥ was significantly different 
fron % , the turbulent rate of change term of equations (3%11) was 
completely negligible, Due to the inability to smooth out the turbulent 
fluctuations of temperature, the definition of the surface layer in this 
study is the layer in which the lapse rate is significantly different from 
dry adiabatic, both the eust and instrumental crrors being considered, Thus 
the surface layer at 1300 LCT extends to about 80 feet, but at 0200 LCT exe 
tends to at least 235 feet. This is evident from the values of Table 10. 
It must be realized that a negative value of this lapse rate is more signifi- 
eimtly different fron vy than a corresponding positive value of the lapse 
rate. Thon, for all practical purposes, equation (3011) within the surface 


laver reduces to: 
oT 2 (oF We 
k GE +x) K, az * 4), Oo, Ge 


where subscrints 1 and 2 refer to a lower and upper level, respectively. 

Let us consider the actual valves of K shovm in Plate I in which the 
value of K is plotted on a-lor-log scale acainst the heicht of the midpoint 
of the slab to which the value of K corresponds. It is noted that, for any 
sounding, the first three points in order of ascending height are very close 
to a straicht line, there being some deviation from the strairht line ro- 
lationship as the fourth, fifth, etc., points are considered, This deviation, 
ecommencins with the fourth noint, is due mainly to two factors: 

(1) Decreasing, therefore less representativo, values of the lapse rate. 

(2) The change in temperature elements above the three foot level; there 


being non-serated resistors below and aerated resistors above. 


(43) 


'«¢ 4 "Te ° 


best line to each 
s 


the maximum height for which K and 2 both in- 
ii 


est fit for each sounding was draym in 


4 


: 
gq 2 . =e: 
the following principles: 
~ ia 


(a) Through the first point. 
(b>) Throush as many additional noints as vossible. 
(c) Separating aprroximately equal numbers of points. 
The slopes and heights to which these lines may be considered representa- 


tive are shown in Table 11 for each soundings 


Vaximum heicht Values of Ky 
Sounding of linear distri- in equation (5.3) 
(Time in LCT) Slope of line, m betion (feet) (em?= sect) 
0200 1.14 2305 365 
0500 a] 2305 506 
1100 Dewi] 9 4540 536 
1300 1659 45.0 2660 
1600 170 2305 410 
1900 1205 4520 152 
2200 1,23 1640 544 


Slopes of Log K Against Log 2 for the Lines of Plate I 


TABLE ll. 
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There is a tendency for the lines for the 1900 and 2200 LCT soundings 
to continue to be representative up to 235 feet, However, for the 0200 and 
0500 LCT soundings beyond 23.5 feet, the points fall well below the line» 
This sugsests the possibility of fitting a curve other than a straight line 
(perhaps a parabola) in the early morning hours; however, this idea was not 


pursued further. 


If we write the equation of any straight line of Plate I in the form 
los K = m(log Z) +d, (502) 


where m is the measured slope (see Table 11) of the line, we obtain upon 


raisins to powers of 10, the equation 
K = 10° 2" K tam (503) 


Since each line passes throuzh the point corresponding to Z = 1/4 footy the 
value of K, for each soundins may be assisned by setting Z = 1/4 and K to 
its appropriate value. The values of I, for cach sounding are shorm in the 
last column of Table 11. Eouation (563), with the appropriate values of Ky 
and m, represent the distribution of k with respect to height for the various 
times of day and un to the heights indicated in Table lle This is the chief 
empirical result of this study. 

Finally, we shall state the conclusions resulting from the forecoing 
analysis end from Pigures 4 throuch 6. These conclusions are: 

(1) The averace mid-day value of m is 1.73. 

(2) The averare value of m near midnirht is 1.18. 

(3) The value of m seems to approach unity at sunset. 


(4) The exponent m appears to be a function of stability. 
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where Z is the heirht in feet, 
(6) The minimm value for K, level for level up to 60 feet, appears 
to occur about 1900 LCT, with a very slow rate of increase during the night. 
(7) From Picure 6 at 125 feet there apvears to be a second pronounced 
dinrnal maximum of K at 2200 LCT. This appears to be due to the diurnal 
variation in the heisht of the surface layer. According to Lettau fs] ; 
the surface layer is characterized by increasins values of kh fith height, 
with decreasing values of K occurring above this layer. If, during the day-= 
light hours the 125 foot level is above the surface layer, but passes into 
it ‘at about 2200 LCT, the large increase in the value of ¥ at 125 feet ap- 


pears to be reasonable. 
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